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[ By the quenched lattice QCD simulation for two nucleons with finite scattering en- 

QQ . ergy, validity of the derivative expansion of the general nucleon-nucleon potential U(r,r') = 

V{r, \/r)S^(r — r') is studied. The relative kinetic energy between two nucleons is introduced 
through the anti-periodic boundary condition in the spatial directions. On a hypercubic lat- 
tice with the lattice spacing a ~ 0.137 fm and the spatial extent La ~ 4.4 fm with the pion 
. mass rriTr — 530 MeV, the local potentials for two different energies (_B ~ MeV and 45 MeV) 

'— H ' are compared and found to be identical within statistical errors, which validates the local 

approximation of U{r, r') up to -B = 45 MeV for the central and tensor potentials. Central 
^ , potentials in the spin-singlet channel for different orbital angular momentums = and 

■ ^ = 2) at _E ~ 45 MeV are also found to be the same within the errors, which also supports 

the local approximation. 

^ ■ §1. Introduction 

The nucleon-nucleon (NN) potentiaP^^ is a fundamental quantity to study 
various properties of atomic nuclei and nuclear matter. Recently, a first attempt 
to calculate the NN potential from QCD was reported on the basis of the Nambu- 
(<"^ I Bethe-Salpeter (NBS) wave function for the two nucleons on the lattice.'^''^ Also, the 

method has been extended to the baryon-baryon (BB) interactions with strangeness 
the three-nucleon interactiorP and meson-baryon interactions.!^ Since the NN 
interaction is short ranged, the NN potential extracted from lattice QCD simulations 
^ ■ is exponentially insensitive to the spatial lattice extent Lg as long as ^ 

, Then one can calculate observables such as the scattering phase shifts by employing 

the lattice NN potential and solving the Schrodinger equation in the infinite volume. 

In general, the lattice NN potential obtained from the NBS wave function is 
energy-independent but non-local, U{r,r'). In practice, U is rewritten in terms of 
an infinite set of energy-independent local potentials V^^'^\r), V^^^'^\r), ■ ■ ■ , by 
the derivative expansion. These local potentials are determined successively by mea- 
suring the NBS wave functions for different scattering energies E below the inelastic 
threshold Ei^. A possible criterion for the validity of the derivative expansion at low 
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energies is the stability of the local potentials against the variation of the scattering 
energy in the interval < £^ < -Eth- E3 

The purpose of this paper is to check such stability through the lattice data 
at ~ MeV and ~ 45 MeV: these two cases are realized on the lattice by 
taking the periodic and anti-periodic boundary conditions in the spatial directions. 
We carry out quenched lattice QCD simulations with Lg ~ 4 fm and the pion mass 

~ 530 MeV. We will show that the leading-order local potentials at the above 
two different energies show no difference within statistical error, which validates 
the local approximation up to = 45 MeV for the central and tensor potentials. 
**)| Difference of the spin-singlet central potentials between £ = and ^ = 2 is also 
studied, with £ being the orbital angular momentum. A preliminary account of these 
results is given in RefsJT3 ]l . [H|) . 

This paper is organized as follows. In Seci2l we make a brief review on the 
energy-independent non-local potential and its derivative expansion. An explicit 
construction of the leading order terms of the derivative expansion is also presented. 
In SeciSl we explain a method to realize non-zero energy NN scattering on the lattice 
through the spatial boundary conditions. In particular, we introduce a novel momen- 
tum wall source operators which are suitable for the purpose of the present paper. In 
SeclU we present numerical results for the NBS wave functions and the associated 
leading order potentials for different E and £. SeciS] is devoted to summary and 
concluding remarks. In Appendix [XJ we give a brief summary of the representation 
of the cubic group used in this paper. In Appendix [Bl some details of constructing 
the 1 = 2 source operator by using the cubic group representation is presented. 

§2. Non-local NN potential and its derivative expansion 

To define the NN potential in QCD, we consider the equal-time Nambu-Bethe- 
Salpeter (NBS) wave function in the center of mass (CM) frame defined by 

^o.p{r-k) = {{)\p^{x)np{y)\B = 2-k), {r = x-y), (2-1) 

where \B = 2;k) is a QCD eigenstate with baryon number two (B = 2), and 
Pa{x), np{y) are local composite nucleon operators with spinor indices a and /3. 
The asymptotic relative momentum k is related to the relativistic total energy W as 

W = 2^J mj^ + fc^ with mjv being the nucleon mass. In the following, we consider 
the elastic region where W < Wth = 2m]sf + is satisfied with the pion mass m^. 

The asymptotic behavior of the NBS wave function for |r| > R {R being the typi- 
cal interaction range) is characterized by the scattering phase shift for hadrons.lSl'''iSJ^l3 
On the other hand, from the NBS wave function for |r| < R, we can define a k- 
dependent local potential Uk{r) and derive an associated /c-independent non-local 
potential U{r,r') with the use of the information of the NBS wave functions for 

*■* Note that Eth is an observable determined by the pion mass and is considerably smaller than 
the scale of the lattice cutoff . 

**^ In Ising field theory, it is analytically shown that the energy-dependence is weak at low energy, 
indicating that the non-locality of the potential is weak.'^ 
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E < Eth- 

{Vl + e) ^{r; k) = 2fiUk{rmr; k) (2-2) 
= 2^y" U{r, r') (Pir'; k), (2-3) 



where /u = mi^ /2 denotes the reduced mass of the NN system. Derivation of Eq. (j2-3p 
from Eq. ()2-2p is given exphcitly in Ref. [S]). Note also that an equivalence theorem 
between Uk{r) and U{r,r') has been proved in a different manner in Ref. [20]) . In 
practical applications, U{r,r') has an advantages over Uk{r); its /c-independence 
leads to the standard eigenvalue problem for the NBS wave function. Furthermore 
its non-locality can be treated by the derivative expansion, U (r, r') = V{r, \7r)S^{f — 
r'), with 

V{ryr) = Vo{r) + V^{r){cTi-cT2) + VT{r)Si2 + VLs{r) L-S+0{V^), (2-4) 

^ V ' ^ V ' 

LO NLO 

where S12 = 3(<Ti-r)(<T2-r)/r^ — cri-cr2, S = (cri + (T2) /2 and L = rx (— iV^) denote 
the tensor operator, the total spin operator and the orbital angular momentum 
operator, respectively.'^''^ 

Since the total wave function has to be anti-symmetric under the exchange of 
two nucleons, possible combinations of the total isospin /, the total spin S and 
the orbital angular momentum I are restricted to four cases, (/, S, ^) = (l,0,even), 
(0,1, even), (1,1, odd) and (0,0, odd). Thus we may omit the isospin / indices in 
Eq. (|2-4p . Note that the spin-singlet states {S = 0) and spin-triplet states {S = 1) do 
not mix with each other, since the isospin / and the parity P = (—1)^ are conserved 
for QCD with degenerate 2-flavors. To specify two-nucleon scattering states, we 
follow the standard notation, "^^^^ij, with J being the total angular momentum. 

2.1. Spin-singlet potentials 

Let us first consider the spin-singlet channel. Since contributions from S12 and 
L ■ S terms are absent in this case, the Schrodinger equation reads 

{Vl + k^) (Pir; k) = 2^ [Voir) - 3K(r) + {V^, Vp2{r)] + VL2{r)L^ + 0{V^)\ 0(r; k). 

(2-5) 

Terms involving 2n derivatives such as [L^)^ and (V^)" give N^"LO potentials. 
(Note that N^^+^LO potentials are absent in the spin-singlet channel.) At the LO 
level, the Schrodinger equation (|2-5p reduces to 

[Vl + k") cPir; k) = 2iiV^^\r)<P{r- k), (2-6) 

with VQ^^^(r) = Vo{r) — 314- (r). Then the LO central potential in the spin-singlet 
channel is given by 

where E = is the effective kinetic energy between two nucleons. The above LO 
truncation works only when the right hand side of Eq. (|2-7p depends weakly on k 
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and £. This will be checked explicitly in SecS] at low energies and at low angular 
momentums through the comparisons, = 0, -E ~ MeV) vs. = 0, ~ 45 MeV) 
and (£ = 0, ~ 45 MeV) vs. = 2, ~ 45 MeV). 

If k and ^ dependence in the spin-singlet channel becomes visible as these values 
increase, it is a sign of the NNLO terms in Eq. (|2-5p . Then the next step is to 
determine NNLO potentials V^^^'^*^^(r), V^^^'^\r) and V^^^'^^r) through the 
NBS wave functions measured with three different combinations of k and £. Such a 
procedure continues to higher orders as k and ^ further increase. A close analogy of 
this process is the renormalization-scale (k) dependence in the perturvative series of 
quantum field theory; the artificial k dependence of scale-independent quantities is 
canceled order by order as we proceed to higher orders. 

2.2. Spin-triplet potentials 

For the spin-triplet channel, the Schrodinger equation reads 

(V^ + k") ct>{r- k) = 2fi [Voir) + K(r) + VT{r)Si2 + VLs{r)L • S + ^(V^)] <P{r- k), 

(2-8) 

At the LO level, it reduces to 

{Vl + k") ^{r; k) = 2/i + ''^(r)Si2] 0(r; k), (2-9) 

where VQ^*^^(r) = Vo(r) -|- Va{r). 

To be specific, we restrict ourselves to the case with = 1"*" NBS wave function 
to which two partial waves contribute, i.e., ^Si (S-wave [I = 0]) and '^Di (D-wave 
= 2] ) . As shown in Ref . [5|) , Eq. (|2-9p consists of two independent equations 

r(t>{r;k) rSi24>{r;k) \ ( V^^\r) - k^ /2^i \ ( r4>{r;k) 



^(2-10) 

where V (Q) is a projection to the £ = {£ = 2) state. The LO potentials, V^^\r) 

and V^'~'\r), are obtained by solving this 2x2 matrix equation algebraically. 

Spatial symmetry group of the hyper-cubic lattice is the cubic transformation 
group S0(3, Z) instead of the rotation group S0(3, M). Here we employ the = 
representation of the S0(3, Z) for the wave function in the spin-triplet channel \*\ 
Since the spin-triplet belongs to the Ti representation, the = state in general 
contains orbital state R which satisfies G RCg^Ti. Table HIl in appendix I A. II gives 
R = Af , E~^, T2 and . Among them we take the projection to the orbital ^4^ 
representation for V as 



'P^t^ir-k)^^ J2 m->];k), (2-11) 



24 

7^eSO{3,Z) 

where the summation is performed over the cubic group S0(3, Z) with 24 elements. 
The orbital Ai representation is expected to be dominated by the S-wave up to 

*-* Here J is used to represent the quantum number of orbital ® spin even for the discrete group 
SO{3, Z), and P is the parity under the spatial reflection. 
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contamination of higher partial waves with £ > 4. We employ Q = \ — V as a, 
projection to non-A^ orbital components composed of , and representa- 
tions. Non-yl^ orbital components are expected to be dominated by the D-wave up 
to contamination of higher partial waves with £ > A. Note that and contain 
the i = 2 component, whereas does not contain the i = 2 component. 

If k and i dependence in the spin-triplet channel becomes visible as these values 
increase, it is a sign of the NLO terms in Eq. ()2-8p . Then the next step is to determine 
NLO potentials through the NBS wave functions measured with several different 
combinations of k and i. 

§3. Finite-energy NN scattering on the lattice 

To extract the NBS wave function on the lattice, we start with the four-point 
nucleon correlation function, 

Ga(5{x-y,t-to;Jpn) = J^(0|T[p«(a; + r, t)n^(y + r, t) J(to)]|0), (3-1) 

~ cl)^^{x-y;k){B = 2;k\jmO)e-^^'-''^\ t - to > 1, (3-2) 

where the summation over r is performed to select the two nucleon system with 
total spatial momentum zero, J {to) denotes a two-nucleon source located at t = 
whose explicit form will be specified below. The relativistic energy and associated 
asymptotic momentum of the "ground" state of the B = 2 system are denoted by 
W and A;, respectively. As for the sink operators, p{x) and n(x), we employ the 
following local composite operators, 

p{x) = eabc {ul{x)C-f5db{x)) Uc{x), n{x) = eabc {ul {x)Cj5db{x)) dc{x), (3-3) 

where a,b,c are color indices. 

The NBS wave function at ~ MeV is generated under the periodic bound- 
ary condition (PBC), which is imposed on the quark operators along the spatial 
directions. With the PBC, the momentum of a single nucleon is discretized as 
ki = 27rnj/Ls with Ui G Z. Hence, the lowest lying state of the two nucleon system 
in the CM frame roughly corresponds to the state where two nucleons are weakly 
interacting with relative momentum of /cj ~ MeV. The effective kinetic energy of 
such a state is E = k'^/ruN ^ MeV. The NBS wave function at ~ 45 MeV is 
generated under the anti-periodic boundary condition (APBC). Since the nucleon 
also obeys the APBC, the spatial momentum of a single nucleon is discretized as 
ki = {2ni + l)7r/Ls with Ui G Z. Hence, the lowest lying state of the two nucleon sys- 
tem in the CM frame roughly corresponds to the state where two nucleons are weakly 
interacting with relative momentum of ki ~ ztyr/Ls. For the lowest lying state with 
Lg ~ 4.4 fm, the spatial momentum of a nucleon amounts to |fc| ~ ^/^tn/L^ ~ 245 
MeV, which corresponds to ~ 45 MeV in our setup with 1.33 GeV. 

As for the source operators of the two nucleon system, we employ 



Japif) = Paif)Mf), 



(3-4) 
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where Pa{f) and Np{f) associated with a source function f{x) are given as 

P{f) = eabc {Uaif)C^5Dlif)) W), 

N{f) = eabc {Ua{f)C-f,Dl{f)) D,{f). (3-5) 
Here the source operators for u and d quarks are given by 

U{f) = ^(^)/(^)' D{f) ^ d{x)f{x). (3-6) 

X X 

An element TZ of the cubic group S0(3,Z) rotates the quark field operator as 

q{x) ^ q{TZ'^x)A{n~^), (3-7) 

where A denotes the 4-component spinor representation of O as A{e'^) = exp (— jaijLo'^^^ 
with fj^j, = I [7^,7;/]. This leads to the transformation property of Japif) as 

JapU) ^ Ja'fs'in-^ o f)A^,^{n-^)A^,^{n-^), (3-8) 

where (7^-^ o f){x) = f{nx). 

To consider J = and 1, it is convenient to introduce a source operator which 
has definite J and M with J = + 5" and M = + to construct NBS wave 
functions in the ^Sq and ^Si Di channels: 

J^'^^'\f)^Y^ E Jc.p{n-' o f) . p^^p^^''=^\ (3-9) 

7eeso(3,z) 

where Pjf^'^''^ denotes the spin projection operator defined as p^^-O.'^^-O) = (o"2)a/3/\/2, 

P^^"^''^"'"^'^^ = (cr20"M)a/3/\/2 with M = ±1,0, where we take only the upper com- 
ponents of the Dirac indices for simplicity. 

For the PBC, we employ a flat wall (f-wall) source, 

/(f-wall)(^) = 1, (3-10) 

which is invariant under the rotation TZ. Then, Eq. (|3-9p reduces to 

j-{J,Af)^y>(f-wall)-j _ (f- wall) ^^^^ J (f-wall) -J _ p(S=J,Sz=M) ^ (3T1) 

which couples dominantly to the ground state [k = (0, 0, 0)7r/Ls) in the PBC. 

For the APBC, we utilize a set of momentum wall sources ylm-waii) _ {/{*)}j^Q_3 
with 

/(0)(^) =cos((+x + y + z)7r/Ls), 
/W(r) = cos((-x + y + z)'n/U), 
/(2)(^) =cos((-x-y + z)7r/L,), 

/(3)(^) = cos((+x -y + z)'k/L,), (3-12) 
where the cosine function is chosen to create positive parity states. 
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The cubic group acts on these functions as permutation, which is characterized 
by the cubic group representation, © T2 ■ By taking the A^^ part, Eq. (j3-9p 
becomes 

-, 3 _ 

j=0 

which couples dominantly to the ground state (fc = (1, 1, l)7r/Ls) in the APBC. 
Since this source operator is not translational invariant, it is practically important 
to perform a summation over r at the sink side in Eq. (j3-ip to pick up zero spatial 
momentum states. Instead of Eq. ()3-12p . one may choose a simpler cosine-type 
function 

f{r) = cos(7rx/Ls) cos(7ry/Ls) cos(7rz/Ls), (3T4) 

which gives a source operator coupled to the ground state {k = (1, 1, l)7r/Ls) in the 
APBC. However, it receives a contamination from the coupling with the first excited 
state (fc = (3, 1, l)7r/Ls). In contrast, the source operator with Eq. ()3-12p has an 
overlap neither with the first excited state nor the second excited state, and receives 
contamination only from the third excited state {k = (3, 3, 3)7r/Ls). Therefore, signal 
for the ground state is better for Eq. (|3T2p than that for Eq. (|3-14p . 

Since Eq. (j3-12p contains T2 component, it can be also used to generate the 

state in the channel, which is employed to study the £ dependence of VQ'^\r). 
A general projection formula for the source operator in the spin-singlet sector given 
in Eq. ([BT]) leads to 

^(J=2,M)(^(m-wall).^+) ^ -^^^e'^^^ P^{f^^^)N p{f^^^) ■ p^^=0'5^=0) (3.I5) 

i=o 

for the T2 representation, where M takes 2 and ±1 (modulo 4). See Appendix [Bland 
IA.2l for more details. In the actual numerical calculation, we take linear combinations 
of Eq. ()3T5p to make them into real basis as 

jj=2,xy _ jJ=2,M=2 

jj=2,yz — J_ / jJ=2,M=-l _|_ jJ=2,M=l\ 

v2 



§4. Numerical results 



4.1. Lattice setup 

Employing the standard plaquette gauge action on a 32^ x 48 lattice at /3 = 
5.7, quenched gauge configurations are generated by the heat-bath algorithm with 
the over-relaxation. We accumulate 4000 configurations separated by 200 sweeps. 
The standard Wilson quark action is used to calculate quark propagators with the 
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Fig. 1. (Left) The NBS wave function for the spin-singlet and the orbital channel at _E ~ 
MeV with the PBC. (Right) The NBS wave function in the same channel but at £ ~ 45 MeV 
with the APBC. Both wave functions are normalized as <f>(r = 0) = 1. 
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Fig. 2. Same NBS wave functions as in Fig. [T]but as a function of ; 



hopping parameter k = 0.1665. The Dirichlet boundary condition in the temporal 
direction is imposed at t — to = ±24. The nucleon four-point correlation functions 
are measured for both t — to > and t — to < to improve the statistics by using the 
time-reversal and charge conjugation symmetries.!^ Either PBC or APBC is taken 
in the spatial direction: In the former case, we use four sources at to = 0, 8, 16, 24 to 
improve the statistics. These calculations are performed on Blue Gene/L at KEK. 

From the rho meson mass in the chiral limit, the lattice spacing is determined 
to be = 1.44(2) GeV (a ~ 0.137 fm), which leads to Ls = 32a ~ 4.4 fm. Our 
K corresponds to the pion mass m-,^ ~ 0.53 GeV and the nucleon mass rriN 1.33 
GeV.123 After examining the stability of the A^A^ potentials against the variation of 
t — to, we chose the wave functions and potentials at t — to = 9 in all the plots shown 
in this paper. 

4.2. The NBS wave functions 

Fig. [TjLeft) and (Right) show three dimensional plots of the NBS wave functions 
(j){x, y,z = 0) for the spin-singlet and the orbital Af channel (~^ So channel) at 
£^ ~ MeV and at ~ 45 MeV, respectively. We observe that they behave rather 
differently: The wave function for the PBC is almost constant at long distances. 
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Fig. 3. (Left) The NBS wave function in the spin-triplet and the orbital channel at £ ~ MeV 
with the PBC. (Right) The same at £ ~ 45 MeV with the APBC. Both wave functions are 
normalized as 0(r = 0) = 1. 



which indicates that the asymptotic momentum is nearly zero. On the other hand, 
the wave function for the APBC decreases continuously to zero at long distances, 
since the wave function in orbital state must vanish on the boundary in APBC. 
This can be seen, for example, by using a vr rotation around the x-axis followed by 
the spatial reflection as 4){x,y,z) = 4>{x,—y,—z) = (j){—x,y,z) = —4>{Ls — x,y,z), 
which leads to (j){Ls/2, y, z) = —4>{L^/2, y, z) = 0. 

In Fig. [21 the same wave functions as Fig. [1] are plotted as a function of r. 
Violation of rotational symmetry due to the square lattice can be seen explicitly 
through the multi-valuedness of the wave function at large r for the APBC. Shown 
in Fig. [3] are the similar comparison of NBS wave functions between the PBC and 
the APBC in the spin-triplet and the orbital channel (— ^Si channel). 

In Fig. m^Upper), we plot the NBS wave functions for the spin-triplet and the 
orbital Ti^ channel (~^Di channel). They are highly multi- valued as functions of r at 
all distances simply due to the angular dependence of the orbital T2 representation. 
To extract the radial part only, we divide the wave functions by Y2^m{0^<i>) assuming 
that the angular dependence is dominated by the 1 = 2 component. The results are 
shown in Fig. I^Lower): almost single- valued radial wave functions are obtained for 
both PBC and APBC cases. 

4.3. LO potentials for different energies 

The leading order potentials are extracted from the NBS wave functions ac- 
cording to Eq. ()2-7p and Eq. ()2-10p for the spin-singlet and spin-triplet channels, 
respectively. In order to obtain LO potentials, we need to determine the value of 
E = }^ l(2[i) either from the large t behavior of the NN correlation function or the 
large r behavior of the NBS wave function^ It turns out that the values of E 
from both determinations are roughly equal to their free values, i.e., £^ ~ MeV for 
the PBC and ii^ ~ 45 MeV for the APBC, within statistical and systematic errors. 
Therefore, we adopt these free values as characteristic E in extracting the LO central 



*' We note here that a new method to obtain the potentials by using the f-dependent Schrodinger 
equation has been also proposed recently.'^ 
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Fig. 4. (Upper- Left) The NBS wave function Re(j!)^^(r) for the spin-triplet and the orbital T2 
channel at £ ~ MeV with the PBC. (Upper-Right) The same NBS wave function but at 
B ~ 45 MeV with the APBC. (Lower-Left) The NBS wave function (ji^i divided by the spherical 
harmonics ^21(^1 </>)■ (Lower-Right) Same as the left figure but at i? ~ 45 MeV with the APBC. 
Normalization of these wave functions is fixed uniquely once the normalization of the S-wave 
part is fixed as given in Fig[2]and FigO 




Fig. 5. The spin-singlet central potential, Vc°(r), (Left) at ~ MeV and (Right) at £ ~ 45 
MeV . 



potentials. Note that the tensor potential is free from the uncertainty of E as can 
be seen from Eq. (|2-10|) . 

In Fig. [5l we plot the spin-singlet central potentials VQ^(x,y,z = 0) obtained 
from the corresponding NBS wave functions (/>(x, y, z = 0) in Fig. [TJ Although the 
wave functions have different spatial structure for different energies, the potentials 
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Fig. 6. (Left) The LO central potential V'q (r) for the spin-singlet and the orbital channel 
as a function of r at i? ~ 45 MeV (red solid circles) and at i? ~ MeV (blue open circles). 
(Right) The LO central potential V^°^{r) for the spin-singlet channel as a function of r at _E ~ 
45 MeV , determined from the orbital representation (red open circles) and from the T2 
representation (cray solid circles). 




0120 1 2 
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Fig. 7. (Left) The LO central potential Vl^°\r) for the spin-triplet and the orbital ^Si Di 
coupled channel as a function of r. (Right) The LO tensor potential V^'^\r) for the spin- 
triplet and the orbital "^Si Di coupled channel as a function of r. Symbols are same as in 
Fig.H 



are independent of E and localized in space. 

To make more precise comparison, v!^^\r) is plotted function of r for 
~ MeV (blue) and at S ~ 45 MeV (red) in Fig. [6)^Left). Similar comparisons 
are also made for V^^^ (r) and V^^^ (r) in Fig. [71 In all these cases, we find no 
-E-dependence within statistical errors. We therefore conclude that the LO potential 
is a good approximation for U{r, r') in the energy range = - 45 MeV. 

It should be kept in mind that we employ a large pion mass jn^r ~ 0.53 GeV, 
which may be one of the reasons for the small energy dependence of the LO potentials. 
It is therefore important to increase E or decrease m-,^ and check the point where 
NLO contributions become significant. 

4.4. LO potentials for different orbital angular momentum 

As mentioned in Seci3l source functions in Eq. (|3-12p for the APBC generate not 
only the orbital but also the orbital T2 components. Combining these sources 
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appropriately, one can construct the NBS wave function for the spin-singlet and the 
orbital T2 channel D2 state). Therefore the LO central potential Vq^^^t) can 
be extracted also from this wave function. 

In Fig. [6|^Right), vl^^\r) obtained from the orbital channel is compared 

with the same potential determined from the orbital channel at ~ 45 MeV. 
Although statistical errors are large, we observe that the two determinations give 
consistent result. Assuming that the orbital and representations are domi- 
nated by ^ = and 1 = 2 waves, respectively, we here conclude that the LO potential 
in the derivative expansion is a good approximation of U {r,r') for £ < 2 in the spin- 
singlet and positive parity channel. 

§5. Summary and conclusion 

We have studied the validity of derivative expansion of the energy-independent 
non-local NN potential U{r,r') = V{r,Vr)S'^{r — r'), defined from the NBS wave 
function on the lattice. For this purpose, we have carried out quenched lattice QCD 
simulations for the NBS wave functions with the lattice spacing 0.14 fm, spatial 
lattice size 4.4 fm and the pion mass ~ 530 MeV. Relative kinetic energies 
between two nucleons were controlled by employing the periodic and anti-periodic 
boundary conditions for the quark fields in the spatial directions. 

The leading-order potentials obtained at different energies (i? ~ MeV and 45 
MeV) show no difference within statistical errors, which validates the local approx- 
imation of the potential up to = 45 MeV for the central and tensor potentials. 
We have also compared the central potentials in the spin-singlet channel for differ- 
ent orbital angular momentum = and £ = 2) at ~ 45 MeV. The result also 
supports the validity of the local approximation of the potential at this energy. 

In the future it is important to apply the analysis in this report to the general 
baryon-baryon potentials in full QCD for lighter quark masses (smaller inelastic 
threshold -Eth) and for smaller lattice spacing a to investigate the convergence of the 
derivative expansion in realistic situations. 
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Appendix A 

Spatial rotation on the lattice 

A.l. Cubic group SO {3, Z) 

In this Appendix, we present a brief summary on the symmetry of two nucleon 
system on the lattice. Further account on the representations of the cubic group can 
be seen in RefsJ26 p . l27p . Note that the NBS wave functions in higher partial waves 
on the lattice are first discussed by LiischerP^ 

A relation of irreducible representations between S0(3,Z) and S0(3,M) is given 
in table [J for i < 6. For two nucleon, the total spin S becomes 1/2 1/2 = 1 0, 
which corresponds to Ti{S = 1) and Ai{S = 0) of the S0(3, Z). Therefore, the total 
representation J for two nucleon system is determined by the product J = Ri R2, 
where Ri = Ai,A2,E, Ti, T2 for the orbital "angular momentum" while R2 = Ai,Ti 
for the total spin. In table [Til the product Ri R2 is decomposed into the direct 
sum of irreducible representations. For example, if the two nucleon state in the spin- 
triplet {R2 = Ti) belongs to the = representation, the orbital representation 
Ri should satisfy G Cg) Ti). From the table [III solutions to this condition are 
given by Ri = Al,E+, r+ and T:^ . 

Table I. Numbers of each representation of S0(3,Z) which appears in the angular momentum I 
representation of SO(3,R). P = (—1)^ represents an eigenvalue under parity transformation. 



i p 


Ai 


A2 


E 


Ti 


Ta 


0(S) + 


1 














1 (P) - 











1 





2 (D) + 








1 





1 


3(F) - 





1 





1 


1 


4 (G) + 


1 





1 


1 


1 


5(H) - 








1 


2 


1 


6(1) + 


1 


1 


1 


1 


2 



Table II. A decomposition for a product of two irreducible representations, R\®R2, into irreducible 
representations in S0(3, Z). Note that Ri ® R2 = R2 ® Ri by definition. 







A2 




E 




Ti 






Ta 






Ai 


A2 




E 




Ti 






Ta 




A2 


A2 


Ai 




E 




Ta 






Ti 




E 


E 


E 


Ai ffi 


A2®E 




Ti ® Ta 






Ti ® Ta 




Ti 


Ti 


T2 




era 


Ai( 


B S ® Ti e 


DTa 


A2( 


BE®Tie 


DTa 


T2 


T2 


Ti 


Ti 


© Ta 


A2 i 


SEfSTid 


DTa 


Ari 


B ® Ti € 


DTa 



A. 2. The cyclic group C4 

Elements of S0(3, Z) which correspond to the rotation around the z-axis form 
a cyclic group C4 consisting of four elements 

C4 = {e,C4,(c4)2,(c4)=^}, (A-1) 



14 



Keiko Murano, Noriyoshi Ishii, Sinya Aoki, Tetsuo Hatsuda 



where e denotes the identity, and C4 denotes the rotation around the z-axis by 90 
degrees. It has four one-dimensional irreducible representations, i.e., A, B, Ei and 
E2. They are related to the irreducible representations of S0(2,M) labeled by M = 
0,2, +1,-1 (modulo 4), respectively. The representation matrices are summarized 
in Table Hill A relation of irreducible representations between C4 and S0(3, Z) is 
given in Table ITVl 

Table III. Representation matrices of irreducible representations of C4. 





e 


C4 






A 


1 


1 


1 


1 


B 


1 


-1 


1 


-1 


El 


1 


i 


-1 


— i 


E2 


1 




-1 


i 



Table IV. Numbers of each representation of C4 which appear in each representation of SO(3, Z). 





A 


B 


El 


E2 




1 











A2 





1 








E 


1 


1 








Ti 


1 





1 


1 


T2 





1 


1 


1 



Appendix B 

A source operator for the state 

A general projection formula for the source operator in the spin singlet sector is 
given by 

^5=0;r,M(^(0))^| ^ I)a(7^)J.,(7^-lo/(o)).pg=0-^-o), (B.l) 

■KgSO(3,Z) 

where F labels the cubic group representations {Ai, A2,Ti,T2, E), dp denotes the 
dimension of the representation F, D'^^^iJZ) denotes the representation matrix of the 
representation F, and M is a label of the irreducible representations of C4, contained 
in the irreducible representation F of S0(3,Z). (See Table ITVl ) This M corresponds 
to the azimuthal quantum number up to modulo 4 due to the cubic symmetry. 

To derive Eq. ()3T5p . we consider the subgroup H in S0(3,Z), which leaves 
f^^\r) invariant 

if = {7eGSO(3,Z)|7eo/(0) = /(0)}. (B-2) 

H is generated by C2, which corresponds to the rotation around m = (1,-1,0) by 
180 degrees, and by C3, which corresponds to the rotation around n = (1, 1, 1) by 120 
degrees. H consists of six elements, i.e., H = {e, C3, (03)^,02,0203,02(03)^}, where e 
denotes the identity. We decompose S0(3,Z) by right cosets of H as S0(3,Z) = 
UcgC4 ^'^^ where He = {hc\h E H} denotes a right coset of H in G, which can be 
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labeled by elements of C4. We use the coset decomposition to arrange the summation 
in Eq. ([HO]) as 

jS=0;r,A/^y.(0)^ 

= I E E D[^:,ihc)j^,iihcr^ o /(o)) . p(-o.--o) 

CGC4 h<^H 

= E f E "K'C-) ■ 7 E "E,(<:) = /'"') ■ pfr''''^'\ (B-3) 

where we used o /(^) = Z^*^). Eq. ()3-15p is arrived at by noting the following 
three facts (i) D^^pj^^{c) is diagonal, which reduces to a phase factor c'^^^'^^'^Sm'm, 

(ii) = {c4)~^of^^\ where C4 denotes the rotation around the z-axis by 90 degrees, 

(iii) t Eh^H D^mmW = 1 for r = ^1 and T^. 
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